Localized versus extended systems in density-functional theory: some lessons from the 

Kohn-Sham exact exchange potential 
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A long-standing puzzle in density-functional theory is the issue of the long-range behavior of 
the Kohn-Sham exchange-correlation potential at metal surfaces. As an important step towards 
its solution, it is proved here, through a rigurouos asymptotic analysis and accurate numerical 
solution of the Optimized- Effective-Potential integral equation, that the Kohn-Sham exact exchange 
potential decays as \n[z)/z far into the vacuum side of an extended semi-infinite jellium. In contrast 
to the situation in localized systems, like atoms, molecules, and slabs, this dominant contribution 
does not arise from the so-called Slater potential. This exact-exchange result provides a strong 
constraint on the suitability of approximate correlation-energy functionals. 
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Introduction In their seminal density-functional- 
theory (DFT) investigation of the electronic structure 
of metal surfaces, Lang and Kohn [l| pointed out that 
far outside the surface the Kohn-Sham (KS) exchange- 
correlation {xc) potential Vxc{z) of DFT should behave 
like the classical image potential — e^/4z, z being the 
distance from the surface. While from the physical view 
point this suggestion results attractive and reasonable, 
forty years later its rigorous proof is still an open ques- 
tion. Two kind of approaches are possible to address 
this difficult problem. One, followed already by some au- 
thors : 2;-4] , is to consider exchange and correlation contri- 
butions to the KS xc potential together, we note, however, 
that as the correlation contribution should always be ap- 
proximated, great care must be taken in approximating 
the corresponding exchange contribution in a compatible 
way. Within this context, it is not surprising that various 
asymptotics have been suggested for the KS exchange- 
only (x-only) potential along this pathway 043] ■ 

A second way to proceed, followed in the present Let- 
ter, is to exploit the fact that since the exchange- energy 
functional is known exactly, the corresponding KS ex- 
change potential Vx{z) can also be known exactly, by 
using the Optimized-Effective-Potential (OEP) method 
of DFT And knowing the exact 14 (z), the analysis 
of the more elusive KS correlation potential Vc{z) may 
be advanced on firmer grounds than previously. We have 
succeeded along this second type of approach, by prov- 
ing rigorously that the asymptotic behavior of the KS 
exchange potential far into the vacuum side of a semi- 
infinite jellium is of the form ln(z)/z. This analytical 
result is supported by a fully self-consistent numerical 
solution of the OEP integral equation, which describes 



accurately the KS exact exchange potential at the bulk , 
interface , and vacuum regions of our semi- infinite system. 

Basic OEP equations for a metal surface The calcu- 
lations presented below focus on the semi-infinite (SI) 
jellium model of a metal surface, where the discrete char- 
acter of the positive ions inside the metal is replaced by 
a uniform distribution of positive charge (the jellium). 
The positive jellium density is given by n+{z) = n 0{—z), 
which describes a sharp jellium (z < 0) - vacuum (z > 0) 
interface at z = 0. The model is invariant under trans- 
lations in the x,y plane (of normalization area A), so 
the single-particle KS orbitals of DFT can be rigorously 
factorized as (ySkn 

fc(r) = e'^rP^f^{z)/^/AlJ, where p and 
k|| are the in-plane coordinate and wave vector, respec- 
tively, k and L refer to the remaining (continuous) quan- 
tum number and the normalization length, both along 
the z direction. The effective one-dimensional KS spin- 
degenerate orbitals Cfc(z) are the solutions of the KS dif- 
ferential equation (atomic units are used throughout) 
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where eu are the KS eigenvalues, and Vks(-z) = Vh(z) + 
Vxc{z). Vh(z) is the Hartree potential, and Vxc(z) := 
6Exc/Sn{z), with Exc ■= Exc[{ek}, {^k}] and 14c(z) be- 
ing the xc energy functional and potential, respectively, 
and n(z) the ground-state electron density. The OEP 
integral equation whose solution provides the KS xc po- 
tential Vxc{z) is compactly given as Q 



{ki - e)n{^)Uz) dk + 



= 0. 



(2) 



2 



Here, 4'/c(z) are the so-called orbital shifts, defined by 

/•oo 

^k{z)^ ik{z')/W^,{z')Gk{z\z) dz' , (3) 
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being the KS Green function, AVj:^{z) = V^dz) 
and = [47r/A{kl - e)Ckiz)]SE,JSUz); u^z. 

are usually referred to as orbital-dependent xc potentials. 
The symbol "P" in Eq. Q denotes the "principal value" , 
and kfp represents the magnitude of the Fermi wave vec- 
tor [7|. The exact Vxc{z) entering Eq. ^ is obtained 
as the solution of this integral equation, which must be 
solved self-consistently together with Eq. ([1]). 

While for formal discussions the integral form of the 
OEP equation [Eq. ([2])] is useful, it is often more conve- 
nient to recast it in the following fully equivalent form, 
after a well-established sequence of transformations [1,01: 
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where V^^^{z) represents the so-called Krieger-Li-Iafrate 
(KLI) contribution 
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and 



V:. 
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with dk = (kp — fc^) dk, primes denoting derivatives with 
respect to the coordinate z and the ground-state electron 
density n{z) being given by the following expression: 



n{z) 



= ^ j\kl-k')\U^)\^ dk 



(8) 



Now we focus on the exchange contribution Vx{z) to 
the KS xc potential of Eq. ([5]), as obtained by replac- 
ing the orbital-dependent xc potentials u'^^iz) entering 
Eqs. ([6|) and ([7]) by their x-only counterparts which 
are known exactly ^Oj. In this case, the first term on the 
rhs of Eq. ^ is easily recognized to be twice the position- 
dependent exchange energy per particle £x(,z), defined 
as the interaction between a given electron at z and its 
exact- exchange hole Noting that 2ex(z) — V^iz), 

Vx{z) being the so-called Slater potential |8j], we write: 



Vx{z) = V^{z) + V,j^{z) + V^^^'\z) 



(9) 



where V.^{z) represents the contribution to the exchange 
potential Vx{z) from the x-only counterpart of the second 
term on the rhs of Eq. ^ . 




[{kl-k^)^k{z)ik{z) + %{z)i'^{z)\ ~ 
27r2n(z) 



FIG. 1: (Color online) OEP self-consistent calculations of 
14"^ (2) (top panel) and V^^'''\z) (lower panel) for the SI ge- 
ometry (solid line) and jellium slabs of various thicknesses d, 
for rs — 2.07. Xf = 2n/kF is the Fermi wavelength. The 
jellium- vacuum interface is at z = 0. 



Numerical results In the case of the SI jellium, we 
have achieved the self-consistent numerical solution of 
the x-only counterparts of Eqs. ([T]) and ([5]). The KS 
equations have been solved by following the procedure 
explained in Ref. [l^ ; the orbital shifts were directly cal- 
culated from its definition in Eq. (jSj, with the KS Green 
function computed using the procedure of Ref. Til. 

The correct asymptotics (at z — > 00) of tlic Slater po- 
tential Vx{z) have been reported for a SI jellium )3t- 14 
and for jellium slabs Pi|, with the result that V^{z) de- 
cays in both cases as —fi/z, but with a coefficient /3 that 
in the case of a SI jellium is electron-density dependent 
while for jellium slabs /3 = 1. Hence, here we focus on 
the remaining contributions: V^{z) and ^^^^"'(z). 

OEP self-consistent calculations of V^{z) and 
yShiitf^^-^ for an electron-density parameter corre- 
sponding to the average density of valence electrons in 
Al (rs = 2.07) are plotted in Fig. [U for the SI geometry 
and for jellium slabs of various thicknesses. We note that 
the bulk limit is correctly reproduced: V^{z — > —00) = 
/cF/27r ~ 0.148Hartree while V^^'^\z -00) 

presents small oscillations around zero, as it should be. 
In the case of jellium slabs, the numerical results were 
obtained using the procedure followed in Ref. 

In the vacuum, both V^[z) and V^^^^^'iz) decay ex- 
ponentially for jellium slabs, as discussed before [6|. In 
the case of a SI jellium, however, the decay of V^{z) 
and Vx^^^^{z) turn out to be qualitatively different. This 
is seen in Fig. [21 where these quantities are plotted for 
rs = 6, together with the Slater potential V^iz) and the 
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FIG. 2: (Color online) OEP self-consistent calculations of 
Vi{z), Vj^{z), and V^'^'^^z) for the SI geometry and = 6. 
The jellium- vacuum interface is at z = 0. The bulk limits for 
V^{z), V^{z), and Vx{z) are correctly reproduced: V^iz — > 
-oo) = ~3kF/2TT ~ -0.153Hartree V^{z -oo) = 
kp ~ O.OSlHartree, and Vx{z — ^ — oo) = —kp/ii — 
—0.102 Hartree. The analytical asymptotes of Eqs. (|12p and 
(|14|l are shown by dotted lines. An enlarged view of the far- 
vacuum region is exhibited in the inset. 



total Vx{z). This figure shows that in the case of an ex- 
tended SI jellium the asymptotics of Vx {z) are dominated 
by V^{z)^ which at large z is positive! 

Analytical asymptotics In order to determine the ac- 
tual asymptotic behavior oiVx{z) in the vacuum region of 
a SI jellium, we first appeal to the asymptotic analysis of 
the KS orbitals S,k{z) entering Eq. (P). One finds [l3 l [l3 | : 



^fe [z ^ <x>) ^ ^kp {z ^ oo) 



and 
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with a — kp/\/2W, W being the work function, and 
Afc = hp — k. This result is perfectly plausible: for 
zkp ^ 1, the only k'a that matter are those close to kp 
and such that zAk ^ 1, so the "window" for the relevant 
fc's decreases linearly with distance. As for the electron 
density, it is interesting to note that it decays with an 
extra power in the denominator that is absent in the 
case of jellium slabs ; this is due to the fact that as 
zkp ^ 1 the factor {kp — fc^) in the integrand of Eq. (jS]) 
becomes effectively small, while in the case of jellium 
slabs this factor is always finite and the electron-density 
decay is purely exponential. This anticipates that quali- 
tative differences might be expected between the asymp- 
totics of localized (slabs) and extended (SI) systems. 

By introducing Eqs. PH)) and PT|) into the expression 
for the Slater potential [the x-only counterpart of the 



first term _on the rhs of Eq. 
result nj 
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([6])], one finds the known 
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Solamatin and Sahni (SS) [lj| then derived the asymp- 
totic structure oiVx{z) from an approximate form of the 
so-called Sham-Schliiter integral equation relating Vx (z) 
to the non-local Hartree-Fock self-energy and concluded 
that the asymptotics of Vx{z) are embodied by half the 
Slater potential, i.e.: Vx{z — )■ 00) = ^V^{z — >■ 00). SS 
supported their result by applying the definition of Vx (z) 
as the functional derivative of the exchange energy (ex- 
pressed in terms of the Slater potential): 
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VAr) = T;VxHr) + ^ / dr' n{r') 



Sn{r) 



(13) 



and then suggesting that the contribution of the second 
term of Eq. ()13p in the vacuum region is zero in the lead- 
ing order of 1/z. It is well known, however, that the 
second term of Eq. (1131) contains another term of the 



form ^V^{r) [l^, leading, therefore, to an expression for 
Vx{z) that contains the full Slater potential V^{z), and 
not one half of it, as noted by Nastos [13i and correctly 
given in Eq. ^ above. Moreover, here we prove that at 
large z the full Vx{z) of Eq. ^ is not dominated by the 
Slater potential V^{z), but by V^{z) instead. 

Here we have succeeded, by introducing Eqs. ()10|) and 
()11|1 into the x-only counterpart of the second term on 
the rhs of Eq. to obtain the following neat (positive!) 
expression for the leading contribution of V^{z) to the 
long-range exchange potential: 



Vx^iz ^ 00) 



'^^ AV'' 2 ~ 

02/" ^Mz^oo\ dk 

g 2Tr^n[z — 00) 

-^[\n{akpz) + C] , 
zTraz 



(14) 



where C ^ 0.96351. In passing from the first to 
the second line we have replaced AVj^'lz), which en- 
ters the calculation of AV^, by its bulk value. That 
is, AV^{z) ~ AV^iz -00) = -kp/ir - ul{z -00). 
The explicit (analytic) expression for u^{z —00) is 
obtained through a k|| Fourier transform of the orbital- 
dependent exchange potential of a three-dimensional 
electron gas [lij . 

Asfor V;S'^'f*(z), we have first analyzed the asymptotics 
of Eq. ^ and then applied to it the operator h'^c,{z). 
Solving the resulting equation for Vx{z 00), one re- 
covers the asymptotic expressions for V^{z) and V^{z) 
given by Eqs. (fT2|) and (fT4|) . and one also obtains: 
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At this point, we need ^'^(^ — >■ oo), which we obtain from 
the asymptotics of the orbital-shifts differential equation 



^fc(z). Noting that 



AK'=(z)-Ay, 

at zkp ^ 1 all contributions in Al/f (z) tend to zero, the 
asymptotics of the orbital shifts are found to be given by 
the following expression: 



-^oo) ^ [f{k) + zg{k)](,k{z oo) 



(16) 



where the first and second terms on the rhs are, re- 
spectively, the homogeneous and particular solutions of 
the orbital-shifts differential equation at z — > oo. Here, 
2g{k) = —AV^/{V2W + aAk), and the explicit expres- 
sion for /(fc) is not needed. After introduction of Eq. 
into Eq. p^ . we find that V^^^^^{z — >• oo) decays as 
ln(z)/z^. Hence, putting this together with Eqs. 
and (jl4l) . we conclude that far outside a semi-infinite jel- 
lium the KS exact exchange potential decays as follows 



K(z -> oo) 



hi{akpz) 
27raz 



(17) 



Equation (|17p represents the main result of this work. 
The asymptotics of Eqs. (fT2|) and (fT4|) are plotted in 
Fig. m (dotted lines), showing that they are in excel- 
lent quantitative agreement with our fully-self-consistent 
OEP numerical calculations at z > Xp . In retrospective, 
the result of Eq. ([TT]) looks natural for the SI case. For 
slabs, Eq. ([9]) yields Vx{z — > oo) = —1/z + AV^ , m being 
the quantum number corresponding to the highest occu- 
pied slab level @ ; the first contribution is brought by the 
slab Slater potential V^{z), and the second contribution 
[brought by the slab V2,^(z)] is a constant which is chosen 
to be zero 0, @1- In the SI case, however, while it is still 
true that the quantity AV^^ entering Eq. is zero, the 
non-negligible contribution from AV^ at kp — k < 1/z 
yields a V^{z) potential that decays as ln(z)/2; and dom- 
inates the asymptotics of the full Vx{z). 

Conclusions In summary, we have solved a long- 
standing problem relative to the long-range behavior of 
the KS exact exchange potential at metal surfaces, as 
an important step towards the understanding of the ac- 
tual asymptotic behavior of the full KS xc potential. 
Through a rigurouos asymptotic analysis and accurate 
numerical solution of the OEP integral equation, we have 
shown that far into the vacuum side of a semi-infinite jel- 
lium the KS exact exchange potential decays as ln(z) / z 
(positive!). This analytical result, which does not arise 
from the Slater potential and is supported by a fully self- 
consistent numerical solution of the OEP integral equa- 
tion, is in contrast to the situation in localized systems, 
like atoms, molecules, and slabs; as in the case of fi- 
nite systems, for jellium slabs the asymptotics of the KS 
exchange potential arise from the full Slater potential, 
which decays as —1/z @. Finally, we note that due to 
the fact that the full KS xc potential of a semi-infinite 



metal should be expected to be absent of the dominant 
ln(2)/z exchange asymptotics, our exact- exchange result 
provides a strong constraint on the suitability of approx- 
imate correlation-energy functionals. 
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